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ABSTRACT 

The shape of the OB-star spectral energy distribution is a critical component in many diagnostics of 
the ISM and galaxy properties. We use single-star HII regions from the LMC to quantitatively examine 
the ionizing SEDs from widely available CoStar, TLUSTY, and WM-basic atmosphere grids. We 
evaluate the stellar atmosphere models by matching the emission-line spectra that they predict from 
CLOUDY photoionization simulations with those observed from the nebulae. The atmosphere models 
are able to reproduce the observed optical nebular line ratios, except at the highest energy transitions 
> 40 eV, assuming that the gas distribution is non-uniform. Overall we find that simulations using 
WM-basic produce the best agreement with the observed line ratios. The rate of ionizing photons 
produced by the model SEDs is consistent with the rate derived from the Ho luminosity for standard, 
log(g) = 4.0 models adopted from the atmosphere grids. However, there is a systematic offset between 
Qo from different atmosphere models that is correlated with the relative hardness of the SEDs. In 
general WM-basic and TLUSTY atmosphere models predict similar T c g, while CoStar predicts T c g 
that are cooler by a few thousand degrees. We compare our effective temperatures, which depend on 
the nebular ionization balance, to conventional photospheric-based calibrations from the literature. 
We suggest that in the future, SpT-T e gf calibrations can be constructed from nebular data. 
Subject headings: galaxies: LMC - HII regions - radiative transfer - stars: atmospheres - stars: early 
type - stars: fundamental parameters - stars: massive 



1. INTRODUCTION 

Through the effects of radiative, mechanical and chem- 
ical feedback, massive stars play a critical role in shaping 
the galaxies in which they reside. These stars ionize their 
local interstellar medium (ISM) and form bright Hll re- 
gions that can be observed in distant galaxies, even where 
the individual stars are not resolved. Since the Hll region 
emission-line spectrum is sensitive to both the spectral 
energy distribution (SED) of the ionizing source and the 
properties of the nearby gas, they are used as diagnostics 
for galaxy pro perties, such as star formation rates and 



wavelengths that are relatively inaccessible to observa- 
tions, we are dependent on the predictions generated 
by stellar atmosphere models to describe the properties 
and SEDs of massive stars. To reproduce the stellar 
SED, the models need to take into account the effects of 
non- LTE conditions, stellar winds, and line-blanketing 
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histories (e.g., Leithercr et al. 1999 Hunter & Massey 



1990 Kennicutt et al.||2000[ ), properties ot the ionizing 



stellar population and the slope of the IMF (e.g.,|Baldwin 
et al.|1981||Kaler|1978t|Rigby fc Rieke|2004| |Stasinska fc 
Leitherer|1996||Copetti et al |1986||Dufour|1975p , and the 
chem ical abundances and chemical evolution of galaxies 
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The shape of the stellar SED is particularly important 
for these diagnostics. For example, it determines the rate 
of ionizing photons (Qo) emitted by massive stars, upon 
which the co mmonly used Ha star formation rate indica- 
tor depends (Kennicutt |1983 ). Other popular diagnos- 
tics use flux ratios of lines with different ionization po- 
tential as diagnostics for th e effective temperature (T e g 



of the ionizing star s (e.g., Vflchez fc Pagel||l988 Stoy 



1933 Zanstra 1927). Thus, the shapes ot massive star 



both challenging and computationally time-consuming 
Therefore, the model atmospheres currently available in- 
corporate these processes with different approximations, 
by balancing an exact treatment of the physics against 
faster computation. The different treatments affect the 
shape of the SED and the properties derived from them 
(e.g., | Sim6n-D faz fc Stasinska||2008| |Voges et aL"1|2008 
Giveon et al. 2002). It is therefore crucial to know how 
well these different atmosphere models represent the true 
SED of these stars. 

Hll regions can be used as a test of these atmosphere 
models. As mentioned earlier, Hll region emission-line 
spectra strongly depend on the shape of the ionizing 
SED. Thus, a direct comparison between the emission 
lines from observed Hll regions and those predicted by 
photoionization simulations will reveal how well the dif- 
ferent atmosphe re models repres e nt the SEDs of mas- 
sive stars ( e.g., |Oey et al .||2000[ |Morisset et al.||2004 



SEDs play key roles in deriving the physical properties 
of nebulae and stars from observed Hn region spectra. 
Since massive stars radiate most of their flux at FUV 
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Esteban et al ||1993| |Crowther et al.||1999| |Stasihska fc 
Schaerer||1997p . For example, previous studies of both 
Galactic and Magellanic Clouds HII regions have shown 
that plane-parallel atmospheres arc too soft and that line 
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this comparison, the metallicity and morphology of the 
gas in the nebula must be constrained. Furthermore, 
most studies use nebulae that contain many stars that 
are distributed throughout the ionized region. Properly 
accounting for these multiple ionizing sources is a ma- 
jor challenge to this approach and vastly limit s the con- 
straints th at one can put on the SEDs (e.g., Ercolano 
et al.|2007 ). Our study circumvents this challenge by us- 
ing single-star Hn regions in the nearby Large Magellanic 
Cloud (LMC). The Hn regions are spatially resolved and 
spherical, which makes the modeling more straightfor- 
ward. We use photoionization simulations to evaluate 
how well atmosphere models represent the shape of the 
massive star SED. 

1.1. Description of the model atmosphere codes 

Our goal in this work is to understand how well the 
OB atmosphere-model grids that are available to general 
users represent ionizing stars. For this reason, we use 
publicly available grids, rather than fitting each stellar 
spectrum in detail to obtain th e appropriate mod el. We 
use the O-star grid presented in Smith et al. ( 2002 ) , here- 
after SNC0 2, in additi on to theCoStar (j Schacrcr fe de| 
Koter|1997|), TLUSTY fLanz fc Hubeny|2003| ) , and WM- 
basic (Pauldra ch" et al.|200lj grids that are already avail- 
able in the stella r atmosphere library of CLOUDY (Fer- 
land et al. 1998). These massive-star atmosphere codes 
include important, but complex, physical processes such 
as non-LTE conditions and the effects of metal lines and 
winds on the transmitted spectrum. Including these pro- 
cesses in detail is computationally expensive. Therefore, 
research groups use various methods and algorithms to 
approximate some of these processes. The differences in 
the algorithms used result in non-negligible differences in 
the SEDs produced by differe nt atmosphere codes (e.g., 
Simon-Diaz fc Stasinska||2008[ ) . 

All the atmosphere codes considered in this work solve 
for non-LTE radiation transfer. Typically this is accom- 
plished by grouping lines of similar excitation energies 
together and applying the same non-LTE correction to 



the populations in tha t group (e.g., Hubeny & Lanz 1995 
Pauldrach et al.pOOl ). In some atmosphere codes, such 
as CoStar and WM-basic, the non-LTE solution includes 
the effects of spherically expanding winds, while others, 
such as TLUSTY, assume a plane-parallel geometry. 

Metal lines in the UV impact the emergent SED in two 
ways, line blocking and line blanketing. Line blocking 
refers to the absorption and scattering of the emergent 
flux due to the higher opacity in the line. Line block- 
ing will increase the temperature in the deeper layers 
of the star because of the scattering of radiation back 
towards the star, a n effect known as backwarming ( [PauT 



drach et al. 20011. Line blanketing refers to the re- 
distribution of energy to r egions where the meta l lines 
are not so densely packed (Pauldrach et al. 2001). The 
main method used to include line blocking and blanket- 
ing is the opacity sampling method. For this method, 
the opacities are evaluated for a grid of frequency points. 
The approximation approaches the exact solution as the 
code increases the number of frequency points it samples. 
CoStar incorporates the opacity sampling using ~ 50 A- 
wide bands in the Monte Ca rlo radiative transfer solution 
(Schaerer & Schmutz 1994). Of the three codes consid- 
ered here, CoStar has the most approximate treatment. 



The TLUSTY atmosphere code includes 180,000 -200,000 
frequency poi nts in the opacity sampling grid (Lanz & 
Hubeny|2003 ). The WM-basic atmosphere code is solved 
m two parts. First, the radiative transfer is solved with 
a fast approximate treatment that samples ~ 1, 000 fre- 
quency points in the Lyman continuum. This fast solu- 
tion is repeated iteratively and generates starting values 
for the final solution. The final solution consists of fewer 
iterations that solve the radiativ e transfer with an exac t 
treatment of the line-blanketing ( Pauldrach et al.[2 001). 

In addition to line-blanketing by metal lines, stellar 
winds chan ge the shape of the ion izing SED at high en 



ergies (e.g., Sellmaier et al. 1996). This is particularly 



important near the 54.4ev ionization potential of Hen 



(Gabler et al. 1989). Both CoStar and WM-basic cal- 



culate the non-LTE solution for expanding atmospheres. 
In the WM-basic atmosphere code, the standard opac- 
ity sampling is modified to take into account line shifts 
due to the expanding winds. These line shifts effectively 
increase t he frequency range tha t can be blocked by a 
given line ( Pauldrach et al.||2001 ). 

While SNC02 is generated with the WM-basic code, 
they use a different set of assumptions for the stellar pa- 
rameters in calculating the SEDs. Specifically, they de 
termine the appropriate mass-loss rate and termi nal ve 
locities a pr i ori using the empirica l rel ations from Prinja 
et al | jl990[ ), |Lamers et al.| ( |1995[ ) and |Kudritzki fc Puis 
( |2000p . To ensure the final atmospheres have these pa- 
rameters, they manually adjust the radiative accelera- 
tion by changing t he force multipliers used in the code 
(Smith et al. 2002). In contrast, the WM-basic grid in- 
cluded in CLOUDY obtains the the mass-loss rate and 
wind terminal vel ocity from the solutio n of the atmo- 
sphere code itself ( Pauldrach et al.||1998 ). We note that 
the SNC02 grid was previously implemented in STAR- 
BURST99 QLeitherer et"allp99 l. 



2. OBSERVATIONS AND METHOD 

2.1. Observations 

To evaluate the atmosphere grids described above, we 
first construct a sample of single-star Hil regions. We 
use the narrow-band emission-line images from th e Mag- 
ellanic Clou ds Emission Line Survey (MCELS; Smith 
et al. 2005) to select small, spherical Hil regions that 
are likely ionized by a single star. Table [ll lists the prop- 
erties of these Hil regions. The first fi ve columns give the 
MCELS ( |Pellegrini et al.|2012[ ), DEM ( IDavies et al.]l976 l 
and stellar designations, and the (J2000) positions. The 
next three columns list observed properties of the ion- 
izing stars. The ninth column has the observed radius 
of the Hn region in parsecs as measured from the Ha 
MCELS image. Figure JT] shows both Ha and three-color 
composites in [Oni], [Sii], and Ha of the Hil regions in 
our sample. 

We used the Inamori-Magellan Areal Camera & Spec- 
trograph (IMACS) on the Magellan Baade Telescope at 
Las Campanas Observatory to obtain both long slit spec- 
tra and Bessell B and V images. Our data were collected 
on the nights of 2008 Jan 29-31 using the f/4 configu- 
ration. In this setup, IMACS has an eight-chip mosaic 
CCD that has a total of 8,000x8,000 15 /im pixels, each 
of which corresponds to 0.11". The seeing was good over 
the observing run, resulting in a final spatial resolution 
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Table 1 

Observed Nebular and Stellar Properties 



MCELS 


Other ID a 


Stellar ID C 


R.A. e 
(J2000) 


Dec e 
(J2000) 


SpT 


vf 

mag 


My 
mag 


-^nebula 

pc 


error 
pc 


lo s( L Ha) 


L 28 


DEM L 08c 


8225 


04:52:11.4 


-66:54:29 


05.5 V 


14.75 


-4.28 


6.7 


1.0 


36.75 


L 32 


N4c b 


8229 


04:52:23.3 


-66:55:16 


B0 V 


15.44 


-3.54 


3.9 


1.0 


35.70 


L 35 




8203 


U4:oz:oo.2 


-00:00:42 


rSU lb 


14.16 


-4.69 


6.5 


1.0 


35.99 


L 43 


DEM L 20 


17251 d 


04:53:30.7 


-67:23:21 


08 V 


14.85 


-3.72 


6.1 


1.0 


36.26 


L 52 


DEM L 26 


19696 d 


04:54:12.0 


-68:21:53 


06.5 V 


14.07 


-4.53 


8.5 


1.3 


36.70 


L 344 


DEM L 276 


50092 


05:40:08.5 


-71:11:02 


07-6.5 V((f)) 


13.98 


-5.03 


9.2 


1.3 


36.82 


L 345 


DEM L 278 


50093 


05:40:09.5 


-71:12:27 


B0.5 lab 


13.15 


-6.30 


5.1 


1.0 


36.27 


L 346 


DEM L 275 


45830 


05:40:11.8 


-69:55:01 


09 V» 


14.69 


-4.16 


4.8 


0.7 


35.97 


L 351 


DEM L 281 


43846 


05:40:43.1 


-70:02:30 


06.5 V 


14.13 


-5.51 


11.3 


1.6 


37.17 


L 390 


DEM L 320 


28307 


05:48:02.0 


-69:53:53 


09 V 


14.46 


-4.78 


5.1 


0.8 


36.45 


L 394 


DEM L 324 


28329 


05:48:43.6 


-69:50:39 


09 V 


14.68 


-4.17 


6.5 


1.0 


35.79 




DEM L 283b 


44979 


05:40:48.1 


-69:43:18 


06.5 V((f)) 


15. 16 h 


-5.19 


5.5 


1.0 


36.30 


a Unless otherwise noted, ID from|Davies et al.| 


{1976J 

















lemze| 1 1956 I i 

c Unless o tnerwise * noTea , ID from OGLE-III | |Udalski et al.|2008| > 
d ID from [Masse"y| l |2002| ) 

e R.A. and Dec correspond to that of the ionizing star. 

f Measurement error on V is 0.05 mag and the syst ematic error is discussed in the text. 
g SpT is inferred but not directly observed, see H3.4| 
h Magnitude from the OGLE-III jUdalski et~aT 2008 ) 

Of - 1.6". 

The long slit observations consist of 3 x 1200s expo- 
sures using a 0.7" slit. We use the 600 1/mm grating, 
which has a spectral resolution R ~ 2730 at Ho. The 
data are binned by 2 and 4 in the dispersion and spa- 
tial directions, respectively, which results in A — 0.76 A 
per pixel. For the first night of the observing run the 
spectral coverage is 3700-5900 A, while for the nights of 
2008 Jan 30-31 the wavelength coverage is 3700 - 6740 A 
due to a different grating tilt. 

We use standard IRAfj^] procedures for the data re- 
duction. The spectra are extracted with the IRAF task 
apextract using separate apertures for the nebula and 
star. The extraction includes a local background sub- 
traction. We extract the nebular spectra in two or more 
apertures. These apertures are selected to exclude any 
stars along the slit, and they otherwise span the entire 
nebula. We sum those spectra to obtain the total emis- 
sion from the nebula along the slit. Figures [2] and [3] show 
the stellar and nebular spectra for the objects in our sam- 
ple. The intensities are scaled to an arbitrary value for 
presentation, and the gaps between the chips are assigned 
a value equal to the continuum level. We flux calibrate 
the spectra using standard stars LTT 3218, LTT 17 88, 
LTT 2415, LTT 2754 and EG 21 ( |Hamuy et aLfl994 |. 

We measure the emission line fluxes from the nebular 
spectra (Figure [3]) with the IRAF task splot assuming 
gaussian line profiles. We use the reddening equation: 



J ( A ) = J o( A ) w -c(fM-nH0)) 



7(H/3) 7 (H/3) 



(1) 



to find the reddening coefficient, c(H^). Here, 
/(A)/7(H/3) and I (\) / I (Hp) are the observed and in- 
trinsic Balme r ratios, respectively, and we use the red- 



dening law of Cardelli et al 
F° r the reddening so 



1989) to determine fx and 
ution we set the ratio of 



3 IRAF is distributed by NOAO, which is operated by AURA, 
Inc., under cooperative agreement with the National Science Foun- 
dation. 



total to selective extinctio n, Ry = 3.45, w hich is ap- 



propriate for the LMC ( |Gordon et al]|2003[ ). The de 
reddened line strengths, relative to H/3, and c(H,3) val- 
ues for each nebula are listed Table [2j In Table [2j we 
also include the H/3 flux that we measure in the slit. 
Note that the slit width is 0.7" relative to the 0.5-2' 
sizes of the nebulae. Therefore, the H/3 flux in Ta- 
ble [2] is not representative of the H/3 flux in the neb- 
ula. We calculate the emi ssion-line measurement error a s 
y/al oat N + a c 2 ont EW/A (jGonzalez-Delgado et al|l994 ) 



Here, <7 C ont is the rms in the continuum near the line, N 
is width of the measured line in pixels, EW is the equiv- 
alent width of the line and A is the dispersion of th e 



spectra in A per pixel ( Gonzalez-Delgado et al. 1994). 
The reddening error and the flux calibration error to- 
gether contribute < 10% error for all our sample, and 
~ 5% for most of our sample. These sources of error are 
combined in quadrature to obtain the flux error listed in 
Table [2 

We assign spectral types (SpT) to the stars from the 
re ctified stellar spectra (F i gure [ 2]) following the criteria 
of Walborn & Fitzpatrick ( 1990J T The spectral type as- 
signments are based on the independent spectral typing 
by four individuals and are accurate to half a spectral 
type. The earliest SpT is 05.5, and the latest is B0. 5. 
Most of the stars are luminosity class V, but our sample 
also contains two B supergiants (Table [l}. 

In addition, we collect single exposures of 20s and 10s 
in filters B and V, respectively. We obtain photometry 
from these images using the IRAF task apphot. Ten of 
our tar gets are in common with stars in the OGLE-III 
survey ( [Udalski et aLp008l ). We compare our measured 
V magnitudes to the OGLE-III magnitudes and find a 
mean difference of -0.07 mag with a standard deviation 
of 0.15 mag. We also have two nebulae, MCE LS L 43 
and L 52, whose stars are in common with the |Massey| 
( [2002] ) UBVR survey of the Magellanic Clouds, and we 
find agreement in V to 0.01 and 0. 09 mag respect ively. 
As a further check, we compare the Massey (2002) and 
OGLE III V magnitudes for ~ 50 stars with V between 
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Table 2 

De-reddened Emission Lines Fluxes™ and Derived Properties 









TOO n 1 b 
LZo 111 


error 


L28 n2 


error 


L 32 nl 


error 


L 32 n2 


error 


L 35 n2 


error 


L 35 n2 


error 








On 


3726 


1.503 


0.115 






0.546 


0.061 




















On 


3726 


3.690 


0.390 


3.446 


0.081 


1.464 


0.211 


2.283 


0.168 


3.705 


0.416 


4.407 










On] 3729 


2.187 


0.159 






0.918 


0.074 




















Nem] 3869 


0.103 


0.009 


0.125 


0.026 
























Oml 4363 


0.009 


0.002 


























Hcl 4471 


0.043 


0.003 


0.043 


0.004 






0.042 


0.018 


0.048 


0.005 










[Om] 5007 


2.015 


0.092 


2.141 


0.027 


< 0.10 


0.009 


0.049 


0.009 


0.068 


0.005 


0.069 


0.008 






Hei 5876 






0.118 


0.002 






0.026 


0.003 


0.045 


0.003 


0.037 


0.006 








Oi] 6300 






0.018 


0.002 










0.021 


0.003 


0.015 


0.004 








Sin 


6312 






0.014 


0.001 
























Nn]6584 






0.265 


0.005 






0.360 


0.014 


0.495 


0.017 


0.482 


0.044 








Sn]6716 






0.221 


0.004 






0.479 


0.020 


0.487 


0.017 


0.470 


0.046 








Sn 6731 














0.333 


0.014 




0.325 


0.033 








c(H/3) 


0.053 


0.058 


0.220 


0.062 






0.150 


0.044 


0.198 




0.182 








log(0/H) 


-3.51 


0.10 


























log(N/0) 






-1.34 








-1.20 
















log(He/H) 


-1.06 


0.01 


-1.06 


0.01 






-1.07 


0.02 














H/3 flux d 


15.22 




23.57 




2.156 




2.741 




3.196 




3.313 










L35 nl c 


error 


L35 nl c 


error 


L43 


error 


L52 


error 


L344 nl 


error 


L344 n2 


error 






On 


3726 


1.835 


0.173 


1.306 


0.136 


1.344 


0.101 


1.493 


0.061 


1.358 


0.121 


1.839 


0.135 








On 


3726 


3.190 




3.293 


0.325 


3.692 


0.197 


3.676 




3.299 




3.046 


0.352 








On 


3729 


2.572 


0.240 


1.884 


0.182 


1.949 


0.124 


2.200 


0.076 


1.941 


0.158 


1.207 


0.108 








Nem] 3869 






0.037 


0.014 


0.040 


0.008 


0.033 


0.006 


0.070 


0.009 


0.129 


0.027 








Oral 4363 














0.016 


0.050 














Hei 4471 


0.008 


0.003 


0.011 


0.004 


0.040 


0.005 


0.028 


0.003 


0.042 


0.004 


0.044 


0.008 






[Om] 5007 


0.034 


0.003 


< 0.03 




1.356 


0.032 


1.377 


0.023 


1.623 


0.075 


1.841 


0.043 






Hei 5876 










0.115 


0.005 


0.110 


0.030 






0.121 


0.005 








Oi] 6300 










0.034 


0.008 


0.014 


0.002 
















Sm]6312 










0.016 


0.008 


0.012 


0.002 






0.013 


0.003 








Nn] 6584 










0.258 


0.011 


0.268 


0.005 






0.280 


0.010 








Sn] 6716 










0.276 


0.010 


0.265 


0.005 






0.222 


0.008 








Sir] 6731 














0.181 


0.004 






0.151 


0.006 








:(HB) 


0.128 




0.133 




0.014 


0.036 


0.032 


0.083 


0.484 


0.056 


0.210 


0.014 






log(0/H) 














> -4.10 
















log(N/0) 










-1.34 




-1.40 








-1.27 








log(He/H) 










-1.09 


0.01 


-1.10 


0.01 






-1.05 


0.02 






H/3 flux d 


3.080 




2.384 




4.644 




13.94 




27.81 




15.51 










L345 nl 


error 


L345 n2 


error 


L346 


error 


L351 


error 


L390 


error 


L394 


error 


DEM L283b 


error 




On 


3726 


0.592 


0.046 


0.592 


0.172 


1.543 


0.134 


1.251 


0.179 


1.520 


0.066 


1.667 


0.287 


1.570 


0.174 




On 


3726 


1.487 




1.144 




3.779 


0.201 


3.271 


0.546 


3.919 


0.452 


4.282 


0.958 


3.859 


0.325 




On 


3729 


0.895 


0.062 


0.552 


0.168 


2.133 


0.174 


2.020 


0.253 


2.399 


0.095 


2.616 


0.375 


2.289 


0.275 




Nem] 3869 










0.033 


0.006 


0.072 


0.035 






0.158 


0.060 


0.039 


0.008 




Oral 4363 














0.005 


0.008 










< 0.002 


0.005 


Hei 4471 






< 0.041 


0.015 


0.031 


0.004 


0.036 


0.008 


0.032 


0.003 






0.031 


0.006 


[Om] 5007 


< 0.027 




< 0.012 




0.905 


0.043 


1.573 


0.078 


0.601 


0.009 


1.452 


0.090 


1.260 


0.064 


Hei 5876 


0.004 


0.002 






0.079 


0.005 


0.090 


0.015 


0.086 


0.003 


0.088 


0.011 


0.079 


0.007 




Oi] 


6300 






0.009 


0.005 










0.010 


0.001 


0.013 


0.007 








Shi 


6312 


















0.012 


0.001 












Nul 6584 






0.295 


0.028 










0.265 


0.008 


0.242 


0.034 








Sn] 


6716 






0.491 


0.043 










0.220 


0.007 


0.343 


0.041 








Sn 


6731 






0.349 


0.031 










0.157 


0.005 


0.231 


0.032 








z(HB) 


0.430 


0.048 


0.389 


0.099 


0.034 


0.295 


0.464 


0.099 


0.302 


0.037 


0.148 


0.063 


0.749 


0.102 


log(0/H) 














-3.4 


0.5 














log(N/0) 


















-1.42 




-1.5 








log(He/H) 










-1.21 


0.01 


-1.17 


0.02 


-1.19 


0.02 


-1.20 


0.02 


-1.12 


0.02 


H/3 flux d 


7.552 




6.041 




9.125 




27.24 




19.00 




2.035 




28.10 





a Flux measurements are relative to H/3. 

b Objects with data taken on different nights are labeled such that nl refers to 29 Jan 2008 and n2 refers to either 30 or 31 Jan 2008. 
c There are two slit positions for MCELS L 32 from each of the two nights we observed it. 

d The H/3 fluxes listed here arc the fluxes within the 0.7" slit in units of 10~ 14 erg s _1 cm~ 2 . They are not representative of the H/3 flux from the 
entire nebula. 
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Figure 1. MCELS images for the objects in our sample. In each image, the left side shows the Ha MCELS image, while the right side is a 
three-color composite of the MCELS bands; Red, blue, and green are Ha, [Om] A5007, [Sn] A6720, respectively. North is up and east is to 
the left. The subfigures show the following objects: (a)MCELS L 28, MCELS L 32, and MCELS L 35, (b) MCELS L 351, (c) DEM L 283b, 
(d) MCELS L 344 and MCELS L 345, (e) MCELS L 43, (f) MCELS L 52, (g) MCELS L 390 and MCELS L 394, (h) MCELS L 346 
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13 and 14 mag. The mean difference (T4iassey — Vogleiii) 
and standard deviation are -0.05 and 0.11 mag, respec- 
tively. From this we see that our values are bracketed 
by the published literature values. We obtain absolute 
V magnitudes, My, after correcting our observed V for 
the measured exti nction and assumi ng an LMC distance 
modulus of 18.48 fWesterhmd]|1997[ ). 



2.2. OB companions 

One important consideration is the possible contam- 
ination from OB companions in our sample. The bi- 
nary fraction for massive st ars is thought to be 40-70% 
(e.g., Sana et al. 2009, 2011) and possibly as high as 90% 



( Kiminki &: Kobumicky||2012|) . Any additional OB stars 
in the nebula will both contribute to the ionizing photon 
budget and affect the shape of the ionizing SED. For OB 
stars, Q n changes by 0. 1 to 0.6 dex for each full step in 
SpT (Smith et al. 2002). Thus, the impact of a compan- 



ion on Qo is maximized for equal mass stars and drops 
rapidly with later SpT. 

The available data put some constraints on possible 
companion stars. For the objects in our sample, the 
nearest resolved stars have magnitudes that are > 2 mag 
fainter than our target stars, thus ruling out contribu- 
tions to Qo from companions at large distances from 
the central star. However, we can still have compan- 
ions on scales smaller than our spatial resolution, 1.6" 
which corresponds to ~ 0.5 pc. In fact, two of our ob- 
jects, MCELS L 346 and MCELS L 390, are confirmed 
eclipsing binaries, and we discuss them in detail in §3.4| 

For the rest of the sample, the observed magnitudes 
will include contributions from both stars, if a compan- 
ion is present. Thus, the luminosities derived from our 
observed magnitudes should be the sum of the luminosi- 
ties from all components and will be representative of 
the system. Furthermore, for a binary or companion, 
one would expect to have spectral type diagnostics rep- 
resenting a mixture of the two spectral types, except in 
the case of an equal mass binary. We note that none 
of our stellar spectra show evidence for composite SpT. 
If any of our objects are equal mass binaries, Q from 
the atmosphere models would be around half th at o f Q 
from L(Ha). We discuss Qo in more detail in £3.3 but 
we note here that we do not find support for equai mass 
binaries from the Qo comparison. Thus, while we cannot 
conclusively rule out binaries from the rest of our sam- 
ple, the data available suggest that any binaries present 
will not significantly affect our results. 

2.3. Method 

To evaluate the stellar atmosphere models, we com- 
pare the Hil regions in our sample to the predictions 
of photoionization simulations. As discussed in the In- 
troduction, the line emission from Hn regions depends 
primarily on the metallicity (Z), the ionization param- 
eter (U), and the SED of the ionizing star. When Z 
and U are constrained, differences between the predic- 
tions from photoionization simulations and the observed 
Hil region spectrum can be directly linked to differences 
in the shape of the model SED and that of the actual 
ionizing star. 

Photoionization calculati ons are performed with ver- 
sion 08.00 of CLOUDY (IFerland et al.||1998b. To set 



up the CLOUDY simulations, we need to match the 
nebular abundance and ionization parameter to the ob- 
served nebulae. In MCELS L 28 and L 351, we detect 



the auroral [Oni 
A4959,5007/[Oin 



A4363 line. We use the ratio [Oni] 
A4363 as input for the IRAF task 
temden to obtain the electron temperature, T e . With 
T e in hand, we derive the oxygen abundance (Table [2| 
using the IRAF task ionic. For this calculation we as- 
sume T e ([Oll]) = T e ([Om]). To set the abun dance of S, 
C, Ne , Ar, Si and Fe, we use the relations from |McGaugh| 
( 1991 1, which relate the elemental abundances to that of 
oxygen. For the rest of our objects, [Om] A4363 falls 
in the gap between CCD chips. In these cases, we adop t 
the mean LMC abundances measured by Garnett ( 1999 ). 
For all our objects, except for MCELS L 346, L 351, and 
DEM L 283b, we calculate the nitroge n abundance from 



the ratio [N il] /[On], as described in Perez-Montero & 
Diaz (2005). Those values are shown in Table [2I We 



note that except for MCELS L 345, the mean difference 
between the calculated log(N/0) and the Garnett (1999) 
values is 0.17 dex with a standard deviation of 0.11 dex. 
This is comparable to the uncerta inty in both the rela- 
tion ( |Perez-Montero fc Diaz 2005J) and in the measured 
abundances above. 
The ionization parameter is defined by 



U 



47rcn e i?| 



(Q ne 2 )s, 



(2) 



where Qo is the rate of H-ionizing photons, c is the speed 
of light, n e is the electron density, Rs is the radius of the 
Stromgren sphere and e is the filling factor of the gas. Es- 
sentially, U describes the ionizing photon density relative 
to the gas density. The ionization parameter will depend 
on both the type of star producing the ionizing radiation 
and the gas distribution around the star. As can be seen 
in Figure [l] the nebulae are Stromgren spheres, which 
means we can assume simple spherical geometry for the 
simulated nebulae. 

The uncertain nebular parameters that control U are 
those that describe the radial density profile of the neb- 
ula: the initial (i?i nnC r) and outer (i? ne buia) radii of the 
cloud, hydrogen density (njj), and e. We use a combi- 
nation of emission-line diagnostics and the MCELS Ha 
images to set these values. To constrain the density dis- 
tribution, we examine the MCELS images and the Ha 
and [Om] A5007 line profiles along the slit. The [Om] 
A5007 and Ha profiles indicate that i?i nn0 r is between 
0.1xi? nebula and 0.5xi? nebu ia- We measure i? ne buia from 
the Ha MCELS images of the nebulae and use those val- 
ues to set the inner and outer radii in the CLOUDY 
simulations. To obtain ns, we assume njj = n e , 
since hydrogen is fully ionized in a typical Hil region. 
The [On] AA3726,3729 and [S11] AA6716,6731 density- 
sensitive doublets both indicate that nearly all our ob- 
jects are in the low-density regime (n e < 100 cm~ 3 ) 
below which these diagnostics are no longer sensitive 
to the density. We also estimate the density using the 
Ha emission measure based on the Ha photometry, and 
find n e < 30 cm -3 , which is consistent with our limit 
from the [Sn] and [On] diagnostics. Based on these con- 
straints, we use nn < 100 enr 3 as an upper limit for 
the densities in all our simulations. Addi tionally, the Hil 
regions in our sample are optically thick \ Pellegrini et al. 
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Figure 2. Stellar spectra from the ionizing stars of the Hn regions. Spectral type changes toward later type from the top to bottom. The 
stellar spectra for MCELS L 52 and MCELS L 394 are boxcar smoothed with a smoothing length of 3 pixels. The flat, noiseless regions 
are the chip gaps, for whic h we assigned a value of unity. MCELS L 346 is an eclipsing binary, and the stellar spectrum shown here is not 
the ionizing star, see j )3. 1| 
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2012| . Therefore, we set n# and e such that the ionizing 
photons are absorbed within the observed nebular size 
(Table [l}. Finally all our simulated nebulae include 5 
km s _1 turbulence. 

In addition to the properties discussed above, dust will 
affect the nebula. In particular, photoelectric heating 
from du st can contribute as m uch as 30% of the total 
heating (van Hoof et al. 2004). Therefore, we include 
both graphite and silicate dust in our simulations. We 
adopt LMC gas distributions and a dust-go-gas ratio, 
Ay/N(ET) = 1.2x l 0~ 22 m ag cm 2 , consistent with Wein- 
gartner fc Prairie] (|200l[), wh ich were introduced into 
CLOUDY in |Fellegrini et ai.| ( |2011[ ). Finally, since our 
observations are based on long silt spectroscopy, we im- 
plement the CLOUDY command slit, which predicts 
the line ratios that would be observed through a slit 
across the simulated nebula. 

We run CLOUDY simulations with the same set of 
nebular parameters for each stellar atmosphere model. 
In addition to specifying which atmosphere grid to use 
for the source, we must also assign the luminosity and 
T e g of the star. We calculate stellar luminosity from My 



using 



log L/L = 0.4 x (Mboj,© - (M v + BC)), (3) 



where M 



bol,© 



is the solar bolometric magnitude and 



BC is the the bolometric correction from Martins et al.l 
( 2005] ), except for the case of MCELS L 345, the B0.5 
lab star, where w e use the bolometric correction from 



Crowther et al.| fl2006| ). The SpT-T eff conversion is 
model-dependent and we leave the T e g as a free parame- 
ter in our simulations. We note that since the BC is tem- 
perature dependent, the input stellar luminosity changes 
with T c g to maintain consistency. For all our dwarf stars, 
we assume log(g) = 4.0. 

We determine the best T e g for each star to be the 
one that produces the correct balance of ionizing flux 
at both high and low energies. We accomplish this task 
by plotting the predicted emission-line flux as a function 
of ionization potential. An example is shown in Fig- 
ure |4j in which we plot the results of six simulations 
that differ in T c g. As expected, as T c g decreases from 
43,000 K to 35,000 K, the predictions for [Oil] and [Nil] 
increase, while the predictions for [Oin] and [Nein] de- 
crease. Thus, the overall slope of the points decreases 
with decreasing T e g. We define the optimal model at 
the turnover T c g, where the slope is flat. For the models 
shown in Figure |4j this would be around T o ff=39,000 K 
(yellow circles). We also prioritize models that match 
the observed Hei lines in our spectrum. The ratio of 
the Hei to Hi recombination lines reflects the level of He 
ionization in the nebula. As long as He remains par- 
tially ionized, this ratio will be primarily dependent on 



Kennicutt et al.||2000[ ). 
en A4686 m any ot our 



the SED of the ionizing source 
We note that we do not detect 
spectra. Therefore, we can set an upper limit on T c g by 
requiring a non-detection of that line in the simulations 
as well 



3. PHOTOIONIZATION MODELS 

3.1. Uniform Density Models 

Figure[l]shows that our HII regions, to first approxima- 
tion, are simple Stromgren spheres. We use the method 



described in |2.3| and the models described in §1.1| to 
first generate a grid of model Hn regions with uniform 
densities. Figures [5] and [6] show the results using this 
prescription for MCELS L 28 and L 43, which are ion- 
ized by an 05.5 V and an 08 V star, respectively. These 
figures are representative of the results from the rest of 
the sample. The left panels of the figures show [Oin] 
A5007 vs. [On] A3727, which reflects variations of U. In 
nebulae with high U there will be more [Oin] relative to 
[Oil], and the points will occupy the top left of the plot. 
As U decreases, [On] becomes stronger relative to [Om], 
and the points will move toward the lower right. The 
right panels show [On] A3727 vs. Hei A5876 . This panel 
illustrates how well [On] is predicted at the appro priate 
T c g, which is traced by Hei ( Kennicutt et al.|20 00). The 
solid and dotted lines show how the line ratios change 
when the -Rmner and T e g are changed, respectively. 

It is immediately apparent that these models do not 
reproduce the observed emission-line ratios. In the sim- 
ulations using CoStar and WM-basic atmosphere mod- 
els, the simulation tracks generally run below and to the 
left of the observed values (Figure [5j left panel). This 
indicates that not enough [On] is produced when [Oin] 
is well predicted and suggests that U is too high. Fur- 
thermore, from the righthand plots, we see that simply 
changing the T e g is not a viable solution; the predicted 
[On] flux does not match the observations for the en- 
tire range in T c g. This is the case for all our objects 
except for MCELS L 43 (Figure M and MCELS L 344 
(not shown). In contrast, the simulated nebulae ionized 
by TLUSTY atmospheres have lower ionization parame- 
ters. The softer TLUSTY SED is able to reproduce the 
observed [Oin] and [Oil] lines in MCELS L 28 (Figure^, 
MCELS L 351, and MCELS L 394. For the rest of the ob- 
jects, the TLUSTY SED is too soft, and the simulations 
have an ionization parameter that is too low (Figure km). 

We rule out the adopted metallicity as the source of the 
discrepancy between the simulations and observations, 
as follows. The metallicity of the gas strongly affects 
the observed line ratios. As the metallicity decreases, 
there are fewer metals to provide cooling, and the neb- 
ula is hotter. Thus, the under-prediction of the low ions 
seen above could indicate that the adopted metallicities 
are too low. To explore how this might affect our re- 
sults, we generate a CLOUDY grid for MCELS L 28 in 
which we change the metallicity from log(0/H)= -3.4 
to log(0/H) = -3.9. The results from this grid for the 
[On] and [Oin] line strengths are shown in Figure [71 The 
measured log(0/H) = -3.51 for MCELS L 28. 

We find that changing the metallicity will not recon- 
cile the discrepancy in the low ions. Although increasing 
the metallicity does decrease the p redicted [Oni]/[On] 
ratio for a given T c g (Figu re 7(a)[ ), it does so by de- 
creasing [Oin] (Figure 7(b)[) . Meanwhile, [On] remains 
un der-pr edicted for the whole range of metallicities (Fig- 



ure 17(c)). Therefore, even though changing the metallic- 
ity might bring the overall [Oin]/[On] closer to the ob- 
served, it cannot account for the discrepancies between 
the simulations and the observations for the individual 
lines. We note that decreasing log(0/H) by 0.1 dex, from 
the observed to the mean LMC value, will change the T c g 
required to reproduce the observed emission-line fluxes 
by 500 K. This is comparable to the 250-500 K uncer- 
tainty in determining T c g using the emission-line ratios. 
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Figure 3. Nebular spectra of the Ho regions in our sample for 3600—4500 A. The spectra are scaled to highlight the weaker emission lines. 
In most of the spectra, the chip gap falls on or slightly red-ward of H7. We assign this region with a value equal to the mean continuum 
level. 
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Figure 4. Ratio of the predicted emission-line flux and the ob- 
served emission-line flux for MCELS L 28. The simulations shown 
above are ionized by WM-basic atmosphere models ranging from 
T cff = 35,000 to 43,000 K. The red horizontal lines are at 
±20% and are representative of the observational variance. The 
best-fitting model is defined as the one with the flattest overall 
slope, with most points lying between the red lines. In this case, 
T eff =39,000 K (yellow circles) has the best fit. 

Since we do not have measured metallicities for most of 
the objects in our sample, this 500 K is included in the 
error on the best fitting T G g, discussed below. 

Another possible cause of the general under-prediction 
of the low ions is the nebular structure. In this case, 
there is not enough dense gas in our simulations at large 
radii to receive a diluted radiation field, which is needed 
to reproduce the emission from lines with low ionization 
potential. This could be because i?i nn or is set too close 
to the star, or because the nebula has a clumpy distri- 
bution of denser gas. An increase of -Rinner increases the 
[On] and [Nil] flux, while it simultaneously suppresses 
the [Oni] and [Nem] flux. At face value, this would move 
the predictions closer to the observations. However, even 
if we vary the inner radius from 0.1 to 0.5 i? n ct>uia, we still 
are not able to simultaneously reproduce the two sets of 
lines (Figure [5|. Furthermore, the inner radius is con- 
strained to be < 0.25 -Rnctmia in all our objects, except 
DEM L 283b, based on the spatially resolved [Om] A5007 
line profiles and the MCELS images. Thus, to within the 
~0.5 pc spatial resolution of our spectral data, these Hn 
regions do not appear to be hollow shells in projection. 
Therefore, neither the metallicity nor i?i nne r can explain 
the discrepancy in these uniform density models. In the 
next section, we explore the effects of a clumpy gas dis- 
tribution. 

3.2. Density Fluctuations 

In Hn regions, non- uniform gas and temperature dis- 
tributions are often invoked to explain the observed prop- 
erties. In fact, direct imaging of Hn regions reveals 
complex morphology, which includes density gradients 
and/or clumps and filaments with high density contrast. 
Non-uniform gas density h as long been invoked to explain 



observed parameters (e.g., Osterbrock & Flather 1959) 



One example is the long-standing problem reconciling 
the discrepant values for T e that result from different di- 
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Figure 5. Predicted line strengths of [Om] vs. [Oil] and [Oil] vs. 
Hel for MCELS L 28 (05.5 V). Line strengths are plotted relative 
to H/J. The solid lines represent models with the same T e g, but 
changing i?i nncr . The effective temperatures shown are T e g= 35, 
37, 39, 41, and 43 kK and increase from the bottom towards the 
top and from left to right, for [Om] vs. [On] and [On] vs. Hel, 
respectively. Note, the CoStar plot only goes up to T c ff=41 kK. 
Dotted lines denote models with the same i?i nncr but changing T c g . 
These lines correspond to 0.10, 0.25, and 0.50 -Rncbulai with the 
radius increasing towards the right in [Om] vs. [On] and towards 
the top in Hel v.s. [Oil]. For comparison, the observed value is 
shown by the black point. 

agnostics. This has been of particular concern because 
the uncertainty in T e translates into significant uncer- 
tainty in the derived abundances. One suggestion is that 
this problem arises from temperature fluctuations in the 
ionized gas, w hich are described by the t 2 parameter ( Pe- 
imbcrt 1967). These fluctuations could arise from small 
clumps of dense gas within the nebulae (e.g., Viegas fc 



|Clegg|1994"{|WilliamsJ1992||Liu fc Danziger|1995| or from 
fluctuations in chemical abundance rather than density 
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Table 3 

Parameters of the best fit simulations. 



MCELS 


L 28 


L 32 


L 35 


L 43 


L 52 


L 344 


L 345 


L 346 


L 351 


L 283b a 


L 390 


L 394 


SpT 


05.5 V 


B0 V 


B0 lb 


08 V 


06.5 V 


06.5-7 V((f)) 


B0.5 lab 


09 V d 


06.5 V 


06.5 V((f)) 


09 V 


09 V 


D. 

■* •'inner 


0.25 


0.25 


0.15 6 ,0.06 c 


0.10 


0.10 


0.10 


0.10 


0.25 


0.10 


0.25 


0.10 


0.10 




4 


4 


9 b(j Q c;nc 
^ .O (O. JU 


4.0 


4.0 


1 


2 6-2.9 


4 o 


4.0 


4.0 


4 


4.0 














CoStar 














-'eft ( K <^) 


oy.o 








ou.u 


J / .u 




ou.u 


OU.t) 




OJ.U 


°,7 n 
o / .u 




0.10 


0.10 


0.25 


0.10 


0.25 


0.10 




0.10 


0.10 


0.10 


0.10 


0.10 


7l e 




7n 
/ u 




ou 




fin 




1 no 


^n 


7^ 


7 F ! 


i on 














TLUSTY 














T e g 


42.5 


33.0 


27.5 


39.0 


38.5 


39.5 


22-24 


37.5 


39.5 


39.0 


35.0 


39.0 


e 


0.10 


0.10 


0.10 


0.15 


0.10 


0.10 


0.15 


0.10 


0.10 


0.10 


0.10 


0.10 


n e 


75 


75 


75 


50 


60 


75 


30 


100 


50 


100 


75 


100 














WM- basic 














T cff (kK) 


41.5 


34.0 


31.0 


39.0 


38.75 


39.0 




38.0 


39.0 


39.0 


37.0 


39.0 


e 


0.10 


0.10 


0.15 


0.10 


0.10 


0.10 




0.10 


0.10 


0.10 


0.10 


0.10 


n e 


75 


75 


25 


50 


60 


60 




100 


50 


75 


75 


100 














SNC02 














T cff (kK) 


41.5 


34.0 


31.0 


39.0 


38.5 


39.0 




37.0 


39.0 


39.0 


37.0 


39.0 


e 


0.10 


0.10 


0.15 


0.10 


0.10 


0.10 




0.10 


0.10 


0.10 


0.10 


0.10 


n e 


75 


75 


25 


50 


(50 


60 




100 


50 


75 


75 


100 



a L 283b refers to the DEM ID 

b TLUSTY simulation 

c CoStar, WM-basic, and SNC02 

d SpT is inferred and not observed, see §3.4| 



dRubi n 1989; Kingdon & Fcrland 1995 Giammanco et al. 
|2004[ |Tsamis & Pequi gnot 2005]r~ 

We therefore explore how changing the density struc- 
ture affects line emission by considering a clumpy 
medium. For our analysis here, we assume that dense 
clumps of gas are uniformly sprea d throughout the neb- 
ula, w ith a vacuum between them ( Osterbrock & Flather 
1959|). In reality, the clumps will be interspersed with 

However, the emission mea- 
Therefore, most of the ob- 



lower density, diffuse gas 
sure is proportional to 
served flux is coming from the dense clumps, and we can 
treat the inter-clump space as though it were a vacuum. 
In the simulations, these clumps are treated in a statisti- 
cal sense; both the volume emissivity of the gas and the 
optical depth along the line of sig ht are decreased by e in 
each zone calculated by the code ( Ferland et al.| [T998). 

In our final simulation grid, T e ff and e are varied, with 
the density chosen as discussed above. A decrease in 
the filling factor of the gas necessitates an increase in 
the density of that gas to maintain the nebular radius 
(Equation [2J . In this scenario, dense gas is present at 
small radii to experience a strong radiation field, as re- 
quired for the lines with high ionization potential. At the 
same time, more dense gas is present at larger radii from 
the star that receives a diluted radiation field to repro- 
duce ions with low ionization potential. Our constraints 
on the sizes of the simulated nebulae, combined with our 
upper limit on the density, effectively provide us with 
lower limits for nn and e. We cannot decrease the den- 
sity past the value that leads to correctly sized nebulae 
when e = 1.0, and we cannot decrease e below the value 
that corresponds to n H — 100 cm~ 3 . We note that the 
filling factor set by the density limit, e ~ 0.10 in most of 
our objects, is consist ent with the IP -1 — 10~ 3 commonly 
used i n the literatur e (Kcnni cutt|1984| Giammanco et al. 
20041 iHunt & Hirashita||2009p. Furthermore, the nebu- 



lae in our sample exhibit substructure (Figure [TJ that is 
consistent with this level of density fluctuations. 
Table [3] shows the parameters of the best fitting model 



for each nebula. These parameters include i?i nn cr and 
log(g) for each nebula and star, as well as the T e g, e, and 
riH that correspond to the best simulation for each atmo- 
sphere model. Figure [8] shows the results of these simu- 
lations. In this Figure, we plot the ratio of the predicted 
emission-line flux to the observed, -FpreZ-Fobs, as a f unc_ 
tion of the ionization potential of the emission line. Each 
model atmosphere corresponds to a different color and 
symbol. The error bars represent the observational er- 
ror, which includes the line measurement, reddening cor- 
rection and flux calibration errors. As discussed in S2.3 



these simulations are selected by finding the combination 
of free parameters (e, T e g) that reproduces the observed 
Hei emission-lines, and yields the flattest slope with most 
points lying within ±20% of unity, denoted by the red 
horizontal lines. The 20% takes into account the obser- 
vational variance, which is estimated from MCELS L 28, 
L 32, L 35 and L 344, for which we have mult iple nights 
of da ta a nd / or multiple slit positions (Figures 
8(j)j and |8(e)[ ). The additive error on our se 
and e for a given i?i nnC r is ~ 500 



ected J c ff 
1000 K and 0.1, re- 
spectively. The error for T e ff includes the contribution 
based on the un cert ainty in our metallicity assumptions, 
as discussed in §3.1[ 

While they are included in Figure [S] for completeness, 
we find that the [Ol] A6300 and [Sn] A6716, 6731 emission 
lines are not good diagnostics for distinguishing between 
the SEDs. These lines are primarily dependent on nebu- 
lar conditions, and less so on the particular atmosphere 
chosen. This can be seen in Figure|9j where we have plot- 
ted the predictions for Hei A5876 against those of [Ol] 
and [Sn] from a simulation grid that varies e and T e g for 
MCELS L 28. Given the same nebular conditions, [Ol] 
and [Sn] are insensitive to different ionizing SEDs. This 
is particularly true for [Sn] , which changes by up to 25% 
for a large range in T c g. A similar behavior is observed 
for [Ol], although we do see significant changes in [Ol] 
at high T c ff. This might be explained by their relevant 
excitation mechanisms. In the case of [Sn], the ioniza- 
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Figure 6. Same as Figure[5]but for MCELS L 43 (08 V). The 
inner radii shown are 0.10, 0~z5, and 0.50 R n cbula- The effective 
temperatures are 35, 37, 39, 40 kK. 

tion potential of this line is 10.60 eV, which is below that 
of H. In this regime, the differences between atmosphere 
models are small, and so the line will be more sensitive to 
changes in the gas density and morphology. [Oi] A6300, 
on the other hand, is produced via charge-exchange and, 
in the conditions of a typical nebula, depends crit ically 
on the density of H+ ( Osterbrock & Ferland 2006 ) . 



3.3. Rate of Ionizing Photons 

One important comparison is between the rate of ion- 
izing photons produced by the observed nebulae, Q Q jj q , 

and that of the best fitting atmosphere model, Qo,Sed- 
To calculate Q we obtain Ho photometry from the 

MCELS, Ha emission-line image (Table fil The rate 
of ionizing photons calculated from this photometry is 
listed in column 2 of Table |U We note that the MCELS 
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Figure 7. Predicted lines strengths of [Oin] A5007 and [On] 
A3727 as a function of metallicity are shown for simulations of 
MCELS L 28 ionized by WM-basic atmospheres with T cff = 39- 
43 kK (thin lines). These are compared to the observed values for 
MCELS L 28 (thick dashed line) for which the observed metallicity 
is log(0/H)=-3.51 
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Figure 8. The ratio of the predicted flux to the observed flux is plotted as a function of the ionization potential of the emission line for 
the simulation that has the closest match to our observations. The different symbols correspond to CLOUDY simulations that are ionized 
by different atmosphere models: CoStar, TLUSTY, WM-basic, and SNC02 are denoted by the red stars, green circles, blue squares, and 
yellow triangles, respectively. The filled data points use the observed flux from the night of 2008 Jan 30 or 31, while the hollow data points 
are from 2008 Jan 29. 
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narrowband images are not conti nuum subtracted. As a 
result, the error on Q Q Hq ~ 20% ( [Pellegrini et al.|2012[ ). 

Columns 3-6 show the ratio of Qo,Sed and the measured 

It is immediately apparent that the Qo,Sed for dwarf 
stars is offset from Q jj a by +5, -12, -12, and -37 % on 
average for TLUSTY,' WM-basic, SNC02, and CoStar, 
respectively. This excludes MCELS L394, for which 
Qo,Sed is 3-5 times higher than <5 jj q (see 93. 41), and 
the two supergiants. From this diagnostic, TLUSTY, 
best represents the ionizing source. However, both WM- 
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Figure 9. [Oi] A6300 and [Sir] A6716 vs Hei A6876 for 
MCELS L 28. Dotted lines are models with constant i?; nner =0.10, 
0.25, and 0.50 fl n ebula an d solid lines have constant T e g= 35, 37, 
39 kK, except WM-basic, which also includes T eff = 40 kK. The 
values of i?i nnor increase going upward and T e ff increases from left 
to right. 



basic and SNC02 also do well, particularly within the 
20% error on the Ha photometry. We note that the off- 
set does depend on atmosphere model, but we find no 
trend with spectral type or Ay- 

One important factor that affects this comparison is 
the assumed gravity. In our simulations, we use log(<?)= 
4.0 for all luminosity class V stars. However, these stars 
may have log(g) in the range from 3.8 to 4.1 (e.g., |Her 



rero et al.|1999| |Massey et ai.|2005| |Martins et at| 2005 ) . 

Decreasing log(g) by just 0.1 dex can increase Qo,sed 
by as much as a factor of 2. This change more than 
compensates for any discrepancy shown in Table |4j and 
supports that Qo.Sed is consistent with Q n a , within 
the uncertainty. 

We note that T e g will affect <5o,sed- Both the bolo- 
metric correction and the shape of the SED play a role in 
determining Q , and both properties depend on T c g . As 
noted in j ]3.2| the error on T e g is between 500 and 1000 
K. If T c ff is too low by 1000 K, Qo.sed would increase 
by 10-15%. However, a change in T e g of 1000 K sig- 
nificantly affects the ionization structure of the nebula, 
which worsens many of the fits shown in Figure [HI Thus, 
uncertainty in T e g does not significantly contribute to 
any offset in Qo,SED/<3 0i H a - 

In addition to the general offsets in Qo,Sed/Q jJq,, we 
also find differences between model grids. In general, 
TLUSTY has the highest Qo,sed/Qq Hce with most val- 
ues near unity, and CoStar has the lowest, with many val- 
ues near 0.6. WM-basic and SNC02 fall between them. 
The trend between atmosphere models reflects the im- 
portant role that the relative hardness of the SED plays 
in the temperature selection process, which in turn af- 
fects Qo. For a given effective temperature, differ ent at- 
mosphere models will produce different Qp (e.g., |Voges| 



et al.||2008| |Sini6n-Dfaz fc Stasinska||2008|) . On top ot 
that, we are sensitive to the shape of the ionizing SED 
because our atmosphere-model selection criteria depend 
on matching the emission-line ratios. This means that 
for a softer atmosphere model (e.g., TLUSTY), we select 
simulations with higher effective temperatures than we 
would for simulations using harder atmosphere models 
(e.g., CoStar). This effect leads to the trend seen above: 
the best-fitting models for TLUSTY have higher Qq than 
those of CoStar. Therefore, while the systematic offset in 
Qo is most likely caused by uncertainty in log(g), it may 
also suggest that the WM-basic, SNC02, and CoStar at- 
mosphere models are generally harder than the observed 
stars. However, we note that this is opposite from wh at 



we find based on the [Nem] line, as we discuss in £4.1 

3.4. Individual Objects 

In this section we briefly comment on the individual ob- 
jects in our sample. Unless otherwise noted, we assume 



the mean LMC abundances fromlGarnettl (11999b. Table 

t lists the Helium abundances, which wecalculate from 
e H ei lines using the relations from |Benjamin et al., 
(2002). As a reference, Table [3] shows the parameters of 
our best fitting models 



MCELS L 28 is ionized by OGLE-III 8225 (Udalski 



et al. 


2008) 


sampl 


e, an 



ic m our 
|. In this 



nebula, we detect the auroral [OmJ A4363 emission line, 
from which we derive log(0/H) = —3.51. The Simula- 
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Table 4 

Qo Comparison 



MCELS 


SpT 




CoStar 


Qo,SED 

TLUSTY 


/Qo Ho 
SNC02 


WM-basic 


L 28 


Uil.O V 


48 94 


U.U-L 


u.yo 


U.U 1 


U. 1 1 


L 32 


BO V 


47.91 


1.01 


1.27 


1.08 


0.84 


L 35 


BO lb 


48.20 


0.95 


2.40 


1.62 


0.95 


L 43 


08 V 


48.47 


0.53 


0.97 


0.79 


0.74 


L 52 


06.5 V 


48.90 


0.39 


0.80 


0.64 


0.62 


L 344 


07-6.5 V((f)) 


49.02 


0.54 


0.99 


0.80 


0.75 


L 345 


B0.5 lab 


48.18 




0.16-0.47 






L 346 


09 Y b 


48.48 


1.06 


1.84 


1.27 


1.36 


L 351 


06.5 V 


49.37 


0.27 


0.55 


0.42 


0.39 


L 390 


09 V 


48.65 


0.66 


1.01 


1.16 


1.60 


L 394 


09 V 


47.99 


2.94 


5.36 


4.37 


4.08 


L 283b a 


06.5 V((f)) 


48.51 


0.70 


1.44 


1.10 


1.02 


a Refers to the DEM catalog number 



b Eclipsing binary for which the SpT of the ionizing source inferred and 
not observed, see 33.41 



tions with the best fit have e=0.10 and jih — 75 cm 



As noted in {3.1 based only on the oxygen lines, the 
TLUSTY simulation shows a good fit using e=1.0 and 
riH = 15 cm~ 3 . However, for that simul ation , both 
[Nem] and [Nil] are under-predicted. Figure 
shows this is not the case in the e = 0.10 models 



a) clearly 
If 



we had selected the e = 1 model as our final best fit, 
the corresponding T e g would be 39,000 K. This is 3,500 
K cooler than the adopted T c ff=42,500 K selected with 
the lower filling factor. This T c g difference highlights the 
importance of using as many emission-line diagnostics as 
available to determine the fits. The best fitting T e g are 
39500, 42500, 41500, and 41500 K for CoStar, TLUSTY, 
WM-basic, and SNC02, respectively. 

MC ELS L 52 is ionized by M2002 19696 (|Massey 
2002| an 06.5 V star (Figures \Uj)\ and [8(b)] ). We find 



excellent agreement for lines with ionization potential 
(IP) below 40 eV, and all the atmosphere models pro- 
duce consistent results. For [Nem], whose IP = 40.96, 
this agreement does not hold. The simulation using 
CoStar, the hardest atmosphere, over-predicts [Nem], 
and TLUSTY under-predicts it. Interestingly, we find 
somewhat different nebular parameters produce the best 
fitting simulations between CoStar and the other atmo- 
sphere models. For CoStar, a slightly higher filling factor 
and lower density, 0.25 and 25 cm -3 , respectively, work 
best. Whereas, e = 0.10 and uh = 60 cm" 3 produce 
better fits with the other atmosphere models. We note 
that this difference is close to our error in determining 
e of 0.10, and, thus, might not be significant. The best 
fitting T cS are 36000, 38500, 38750, and 38500 K for 
CoStar, TLUSTY, WM-basic, and SNC02, respectively. 

MCELS L 351 is ionized by OGLE-III 43846, which 
is also an 06.5 V star and, with a radi us of 11 pc, it is the 
largest nebula in our sample (Figures [T(b)1 and[8(c)T). We 
detect [Om] A4363 in this object, from which we calculate 
T e and log(0/H) (Table |2|. However, our [Om] A4363 
detection is at a very poor signal-to-noise, and our result 
is consistent with the mean LMC metallicity. Therefore, 
we use the mean LMC metallicity in the CLOUDY simu- 
lations. As for MCELS L 52, we find excellent agreement 
between simulations and observations for emission lines 
with IP < 40 eV. The best fitting T eff are 36500, 39500, 
39000, and 39000 K for CoStar, TLUSTY, WM-basic, 



and SNC02, respectively. 

DEM L 283b is ioniz ed by OG LE-III 44989, an 06.5 
V((f)) star (Figures 1(c) and [8(d) ). Since we do not have 
IMACS photometry for this object, we use the OGLE- 
III U-band magnitude to calculate the stellar luminosity. 
From the [Om] A5007 line profile, we assign i?i nnor = 0.25 
-Rnebuia- Overall, we find good fits for all atmosphere 
models for lines with IP < 40 eV and e = 0.10. The best 
fitting T cS are 36500, 39000, 39000, and 39000 K for 
CoStar, TLUSTY, WM-basic, and SNC02, respectively. 

MCELS L 344 is ioni zed b y OGLE -III 45830, an 
06.5-7 V((f)) star (Figures [T(d)] and [8(e)] ) . We note that 
the nebular structu re div erges from a Stromgren sphere 
to the east in Figure [l(d) However, our slit position runs 
N-S and covers the regular portion of the nebula. Overall 
we find excellent fits for models where IP < 40 eV. The 
measured [Nem] flux shows a 45% difference between the 
spectra from 2008 Jan 29 a nd 31 , which results in a large 
scatter for [Nem] in Figure 8(e) This apparent discrep- 
ancy is caused by the large observational uncertainty. 
Interestingly, for Jan 31, the simulation using CoStar is 
the only one to match the observations, and the others 
under-predict [Nem] by > 50%. The best fitting T eS are 
37000, 39500, 39000, and 39000 K for CoStar, TLUSTY, 
WM-basic, and SNC02, respectively. 

MCELS L 43 is i onize d by M2002 17251, an 08 V 
star (Figures [1(e)] and |8(f)[ ). As noted in SjO MCELS L 
43 is one of the few objects for which [0inj/]0ii] is well 
predicted by simulations using e = 1 with atmosphere 
models from WMbasic, SNC02 and CoStar. However, 
the right panel of Figure [6] shows that these models can- 
not simultaneously match the Hci lines, which are a more 
reliable T e g indicator. At the temperatures required to 
match Hci, [Om] A5007 is over predicted by a factor of 2- 
3. The fit improves in all emissio n lines when we decrease 
the filling factor to 0.10 (Figure [8(f)] ) ■ The best fitting 
T cff are 37000, 39000, 39000, ancTSSOOO K for CoStar, 
TLUSTY, WM-basic, and SNC02, respectively. 

MCELS L 3 90 is ion ized b y OGLE-III 28307, an 09 
V star (Figures 1(g) and 8(g)). For this object, the sim- 
ulation using CoStar reproduces the relative amounts of 
[Om] and [Nem], while the others do not. This is in 
contrast to the rest of our sample, for which the CoStar 
models are too hard. This star is an eclipsing binary, 
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according to the OGL E-III catalog of eclipsing binaries 



(Graczyk et al. 2011). However, the second star onl 
contributes 10% of the I band flux. Futhermore, Table 
shows agreement between Qo.Sed and Q Therefore 
the second star is not significant. The best fitting T e g are 
35000, 35000, 37000, and 37000 K for CoStar, TLUSTY, 
WM-basic, and SNC02, respectively. 

MCELS L 394 is ionized by OGLE-III 28239, a star 
whose spectral type is between an 09 an d a B1 .5 V, base d 
on the noisy stellar spectrum (Figure 1(g) and |8(h)[ ). 
However, we assign this star the earlier SpT m this range 
because the observed [Nein] rules out stars with SpT of 
B0 V or later. Even with e = 0.10, it is still challenging to 
repro duce the [Sn], [On] and [Nil] lines. As seen in Figure 
8(h) the predictions lie at the edge of our 20% tolerance 
limit. Additionally, all the atmosphere models, except 
for CoStar, do not produce enough flux at the energies 
needed to re produce the observed [Nem] emission-line 
(Figure [8(h)] ). The best fitting T eff are 37000, 39000, 
39000, an7f~39000 K for CoStar, TLUSTY, WM-basic, 
and SNC02, respectively. 

Unlike the other objects, Q Q Ha in MCELS L 394 is a 

factor of a few less than Qo.Sed (Table 5b. This could 
indicate that the nebula is optically thin. If the nebula 
were optically thin, we expect to observe less [Oil] and 
[Sn] than if the nebula were optically thick. Based on 
our modeling, however, we see the opposite. We observe 
too much [On] and [Sn] relative to the optically thick 
model. Therefore, while the Ha photometry may suggest 
an optically thin nebula, it is clear that this does not 
adequately explain all the observations. 

MCELS L 32 is ionized by OGLE-III 8229, a B0 V 
star (Figure 1(a) and 8(i) I. [Nein] is not detected in this 
nebula, and thus provides a firm upper limit for the T e g 
of the ionizing star. Both of the grids generated with 
the WM-basic code include atmosphere models down to 
T c fi= 30,000 K, which is adequate for this B0 star. How- 
ever, the CoStar grid in CLOUDY does not have dwarf B 
star models. The coolest CoStar dwarf, T o gf=34,000 K, is 
still a bit too hard to match the observed emission lines. 
We optimized the T e s to reproduce the Hei lines, with 
the result that [Oin] is under-predicted. The best fitting 
T cS are 34000, 33000, 34000, and 34000 K for CoStar, 
TLUSTY, WM-basic, and SNC02, respectively. 
MCELS L 35 is ionized by OGLE-I II 820 3, one of 



our two supergiants, a B0 lb star (Figure 1(a) and 8(j) ). 
In addition to varying the T e g and e, we vary the surface 
gravity bet ween log(g) = 2.5 — 3 . 5, as appropriate fo r 



supergiants (Crowther et al. 2006 Trundle et al. 2004) 



TLUSTY is the only atmosphere code which includes this 
whole range of gravities and temperatures. The lowest 
gravity in the available grids for WM-basic and SNC02 
are 3.0 and 2.95, respectively. In the case of CoStar, 
the lowest available gravity for T c ff>26,000 K is log(<?) 
= 2.86. For TLUSTY, we find that simulations using 
log(g) = 2.6 — 3.0 all yield good fits, although at dif- 
ferent temperatures. For example, for log(g) = 2.6, the 
best fitting T eff = 27,250 K, whereas for log(g) = 3.0, 
T cS = 28,500 K. For WM-basic, SNC02 and CoStar we 
use log(<?) = 3.5 and find T c g= 31,000 K for all three 
atmospheres. 

MCELS L 345 is ionized by OGLE-III 50093, the 
other supergiant in our sample, a B0. 5 lab star (Figure 



l(d)| and [8(k)[ ). Based on the [Nil] /[Oil] ratio, this object 
has a much higher nitrogen abundance than the rest of 
our sample, log(N/0) = —0.77 compared to the mean 
LMC value of log(N/0) = —1.5. This higher abundance 
is confirmed by the simulations because those that use 
the mean value are unable to reproduce the [Nil] /[On] 
line ratio. [Oin] A5007 is not detected in this object, 
which places an upper limit of T c ff=26,000 K. As dis- 
cussed for MCELS L 32, we are only able to test the 
TLUSTY atmosphere models for this object. Note, the 
points along the bottom of Figure 8(k) are lines for which 
we only have upper limits on the line strength, and the 
simulations are consistent with the data. We find the 
best fitting T eB falls in the range of 22,000 - 24,000 K. 

MCELS L 346 is ionized by an ec lipsing binary, MA- 
CHO 81.9725.16 ( |Alcock et aL|T997| . The spectral type 
determined from our observed data is Bl-1.5 V (Figure 

t. However, we do detect [Nein] in this object, despite 
e late SpT. The strengths of the [Nem] and [Oin] lines 
suggest that this nebula is ionized by a star of much ear- 
lier spectral type. Based on the [OmJ/H/3 and [Nein]/H/3 
line ratios compared to the rest of our sample, the ion- 
izing star in this nebula likely has a SpT between 08 V 
and 09.5 V. This analysis is further supported by the 
measured log(Q Hq )=48.48, which falls into the [47.99- 

48.65] range spanned by the 08 - 9 V stars in our sample. 
Considering the short period of this binary, 1.2 days, it 
is possible that we observed the nebula when the earlier 
type star was eclipsed by a later type companion. Since 
the early type star will dominate the ionizing SED, we 
model this object with a single atmosphere model. 

We note that one side of this nebula has an irregu- 
lar mor pholo gy, while the other is more Stromgren-like 
(Figure | l(h) ). Therefore, we only use the emission-line 



ratios from the portion of the slit that lies along the 
regular half of t he ne bula. The best fitting models are 
shown in Figure |8(1)[ These models have Ri nnC i— 0.25 
-Rnebuia, e = 0.10, and hh = 100 cm~ 3 . We find ex- 
cellent fits for WM-basic, TLUSTY, and SNC02 except 
at [Nem] A3869, while CoStar over-predicts both [Om] 
and [Nem]. The best fitting T off are 36,000, 37,500, 
38,000, and 37,000 K for CoStar, TLUSTY, WM-basic, 
and SNC02, respectively. 

4. DISCUSSION 

4.1. Comparison of Different Atmosphere Models 

In the previous section we discussed the simulations 
for each object individually. However, the question at 
hand is: how well do the atmosphere models reproduce 
the ionizing population, in general? In this section, we 
discuss the trends we find across our sample. To look at 
the results collectively, we plot Fp re /F b s for individual 
emission lines as a function of the spectral type of the 
ionizing star. These plots are shown in Figure |10| 



With the exception of [Nem] A3869 (Figure [l0[d)|), Hn 



region simulations using CoStar are consistent with the 
other atmosphere models and the observations (Figure 



10). For high io nizatio n potential lines, such as [Nem] 



shown in Figure 10(d) simulations with the CoStar at- 
mospheres over-produce the line emission. Other studies 
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a small range of T c r reproduces the observe d n ebular 
emission for most of our range in SpT (Figure 12). This 
implies that the SEDs generated by CoStar areless sen- 
sitive to changes in the atmospheres between stars of 
different SpT. 

As with CoStar, the TLUSTY models reproduce most 
of the emission lines up to the ionization potential of 
Ne+ 2 . At the IP of Ne+ 2 , 41 eV, the modeled Hn region 
spectrum under-predicts the line emission, and it does 
so to a greater degree than the other atmosphere mod- 
els. This suggests that the TLUSTY atmosphere is not 
producing enough flux at high energies and reflects differ- 
encesin the incident SEDs, as shown in Figure [TT] Fig- 



11 compares the SEDs produced by different stellar 
atmosphere codes for two temperatures, T e ff=35,000 K 
and T e ff=41,000 K. At both T e g, the high energy slope of 
the TLUSTY SEP is steeper than the other atmos phere 
models. Similarly, [Simon-Diaz fc Stasihs ka (2008) show 
that for ionizing SEDs with 30, 000 < T eS <~207O00 K, 
TLUSTY produces less flux at high energies than model 
atmo spheres with expandin g winds. Using IR emission 
lines, 
mode 



Morisset et al. (2004) also find that plane-parallel 
atmospheres, 1 LUSTY included, are too soft be- 
tween 27 eV and 41 eV to reproduce the observed nebu- 
lar emission. This confirms that the treatment of the 
expanding wind plays an important role in reproduc- 
ing the flux of ion izing photons with higher energies. 
|Gabler et al.| ( |1989[ ) showed that self-consistent treat- 
ment ot the wind was particularly important for the con- 
tinuum near t he Hell ionizing edge at 54.4 eV. It was 
also shown by Sellmaier et al. ( 1996 ) that the [Nem] 
problem, in which photoionization models consistently 
unde r-predicted [Neiii] A3869 (e.g., Simpson et al.||1995 
1986||Mathis et al.|1985)|Rubin et al.|199ip , is not an is 



sue m simulations using non-LTE atmospheres that have 
a wind extension. We note that while there is a gen- 
eral trend for [Nem] to be under predicted, the scatter is 
large and the degree to which [Nem] is under predicted 
in our TLUSTY simulations is much less than that seen 
for earlier generations of stellar atmosphere models. 

At the opposite energy regime, simulations using 
TLUSTY yield more [On] emission than the other atmo- 
sphere models, and ty pically over- predict [On] 3 72 7 A by 
15-30% (Figure[H)(al While TLUSTY is softer flSimon- 



Diaz & Stasihska 2008b, it is surprising that this trend is 



not seen in the other lines with simila r ionization poten- 
tial, such as [Nil] A6584 (Figure [iO(b)[ ). One explanation 
for this could be the relative shapes of the SEDs above 
the ionization potentials of N+ 2 (29.6 eV) and 0+ 2 (35.1 
eV) , since the flux in that energy range controls the ion- 
ization from 0+ to 0+ 2 and N+ to N+ 2 . For energies 
greater than 35 eV, the luminosity of the TLUSTY SED 
dro ps relative to other models, at least for T e g < 40, 000 
K (Simon-Diaz & Stasihska 2008). Thus, the slight ex- 
cess of [OnJ may result from the softer SED. 

Simulations using atmospheres generated with the 
WM-basic code are the closest to the observed values 
for mo st emission lines (Figure [lOp -c). However, Figure 
|10(d)| shows that the predictions for [Nem] A3869 have 
more scatter than what is seen for the emission lines of 
lower IP. The ratio of predicted and observed flux for 
[Nem] ranges from 0.7 to 1.7, with most points under- 
predicted by 20-30%. As expected, the atmosphere mod- 



els presented in the SNC02 grid are consistent with the 
WM-basic grid implemented in CLOUDY, although we 
find that they occasionally require a slightly different 
T c ff. This is likely due to the different stellar properties 
(stellar radius, mass loss rate and wind terminal veloc- 
ity) used in generating the atmosphere model grid, as 
discussed in fLl] We note that the temperature differ- 
ences, when they occur, are in the range of 250-750 K. 
Although the predictions from the two grids are close, 
the WM-basic grid does slightly better than the SNC02 
grid in reproducing [Om] for late SpT stars. 

4.2. Comparison between SpT and T e g 

The calibration between spectral type and effective 
temperature is an important consequence of the atmo- 
sphere modeling. Ordinarily, the calibration is based 
on careful fitting of photosphcric absorption lines in 
model SEDs to reproduce high-resolution stellar spectra. 
Therefore, the SpT-T e ff calibration depends on the stel- 
lar atmosphere models used in the fitting. This bias is 
especially apparent when one compares the calibrations 
based on atmosph eres codes without line blanketing (e.g., 
Vacca et al. 1996 ) to calib rations based on atmospheres 
with line blanketing (e.g., |Martins et al.|[2005| [Massey 
et al.|[2005 ). The latter calibrations assign effective tem- 
peratures that are several thousand degrees cooler than 
the previous calibrations. 

In contrast to this method, the effective temperatures 
from Table [3] depend on reproducing the ionization bal- 
ance in the nebula that the star ionizes. The points in 
Figure [12] show the effective temperatures that we gener- 
ate as ahmction of the SpTs from Table [l] For compar- 
ison, the lines show calibrations from the literature for 
dwarf stars. The t ypical error on the l iterature calibra- 
tions is ~ 1000 K (j Mart in.-, et~aLl|2005l). Figure [12] does 
not show MCELS L 346 since the SpT of the ionizing 
star is uncertain. MCELS L 394, for which the stellar 
spectrum is noisy, is placed at SpT = 09.5 with a large 
error bar. Based on the nebular spectrum, however, this 
object is probably at the upper end of the shown error 
bar. 

We find that, in general, WM-basic, SNC02, and 
TLUSTY predict similar T e g, while CoStar predicts 
coole r T e ff. This trend is consistent with our findings in 
93.31 the CoStar SEDs are harder than the other atmo- 



sphere models, and therefore require a lower T e g to repro- 
duce the observed emission-line spectra. For the earliest 
stars, the best fitting T c g selected for TLUSTY is slightly 
hotter than that of WM-basic and SNC02. At late SpT 
we see a trend toward the opposite. Since the treat- 
ment of the winds is the biggest difference between the 
TLUSTY and WM-basic atmospheres, this may reflect 
changes in the wind properties toward late SpT stars. 

We note that the two supergiants in our sample have 
T e ff between 24 kK and 31 kK. This large range spans, 
and is consistent with, the expected T e g for Magellanic 
class I star s 25-32 kK (e.g., |Mokiem et al.|2007} |Trundle 
eT^[2004| . 

The effective temperatures presented in Figure [12] arc 
consistent with the relations in the literature, within the 
errors. However, we note that the points appear to follow 
a shallower slope. This is primarily caused by the stars 
of later spectral typ e. The 05.5 - 07 V stars follow the 



LMC calibration of Massey et al. ( 2005 1 (thick-dashed 
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Figure 10. F pre / F obs as a function of spectral type for [Oil] A3727, [No] A6584, [Om] A5007 and [Ncm] A3869 in panels a,b,c and d, 
respectively. Colors and symbols are as for Figure [8] 
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line), but the lat er spectral types pr efer the hotter tem- 
peratures of the Vacca et al. ( 1996 1 calibration (dotted 
line) . 

One contributing factor to the T c g trend is the assumed 
\og(g). We calculate an expected variation in log(g) for 
dwarf stars across our spectral range of 0.05 dex, based 
on the modeled stellar parameters from |Schaerer et al.| 
(1993). If there is an error in our spectroscopically de- 



termmed luminosity class, it will affect the best fitting 
T e g selected. Specifically, if any of the stars classified 
as dwarfs are actually giants or supergiants, the T e g 
will be 1000-3000 K cooler than the T e g chosen for a 
dwarf star. The difference between the log(g) appropri- 
ate for dwarfs and supergiants increases for later SpT. 
Therefore, changes to the SED caused by uncertainty in 
logCq ) will be strongest for stars of late spectral type 
(e.g., Martins et al.|[2005] ). We find 0-1000 K decrease in 
T eff for MCELS L 28, an 05.5 V, going from log(. 9 )=4.0 
to log(fif)=3.0, but a 2000-3000 K decrease in T eS for 
MCELS L 20, an 08 V. 

Another consideration is the mctallicity for which these 
SpT-T ff calibrations are valid, since lower metallicity 
stars will ha ve higher T e g ( |Massey et al. 2005] |Mokiem| 



et al.||2004[ ). All of the literature calibrations, except 
for the thick-dashed line, are calibrations that use solar 
metallicity stars. Thus, we expect most of these cali- 
brations to lie at cooler T e g than the points in Figure 
IT2I While this is the case for the calibrations based on 
line blanketed models, which we see by comparing the 
thick-dashed line to the thin-dashed and solid lines, we 
note the scatter between the points is larger than the 
difference in T e ff between the LMC and the MW calibra- 
tions. Furthermore, the scatter is comparable to the dif- 
ference between calibrations based on atmospheres with 
and without line blanketing. 

The conventional approach to obtaining the SpT-T e fj 
calibration is based on photospheric lines and is sensitive 
to the SED at the photosphere. Our method, however, 
depends on the ionization balance of the nebula and is 
also sensitive to changes in the SED due to layers out- 
side the photosphere. Thus, SpT-T e ff relations derived 
from nebular data have the potential to probe factors in 
the upper layers of the atmosphere models that do not 
match the conditions in the ionizing stars. Our sample 
is small, and as seen earlier, there are a variety of uncer- 
tainties, limiting its usefulness in this regard, but with 
larger samples and better data, this approach promises 
an alternative and important way to refine SpT-T e fj cal- 
ibrations. 

5. CONCLUSIONS 

Diagnostics from Hn region emission-line spectra are 
used extensively to determine the physical conditions of 
the ISM of galaxies, as well as understand their stellar 
populations and their chemical evolution. The shape of 
the ionizing SED is one of the primary properties that 
determines the Hn region emission-line spectrum. How- 
ever, the intrinsic shapes of massive star SEDs are uncer- 
tain, and we depend on atmosphere models to describe 
them. It is therefore crucial to understand and quantify 
how well massive star atmosphere models represent the 
ionizing sources. 

In this study, we compare a sample of single-star 
Hn regions from the nearby LMC to photoionization 



simulations ionized by atmosphere models from widely- 
available CoStar, TLUSTY, and WM-basic atmosphere 
grids. We select our sample of small, Stromgren sphere- 
like Hn regions from narrowband MCELS images. We 
obtain spectra and photometry of this sample using 
IMACS on the Magellan Baade Telescope at Las Cam- 
panas Observatory. We assign spectral types to the ion- 
izing stars in our sample and find a range from 05.5 V 
to B0 V with 10 of our 12 Hn regions ionized by dwarf 
stars. While we choose Hn regions that are likely to be 
ionized by a single star, contamination by OB compan- 
ions cannot be entirely ruled out. Two of the objects 
in our sample, MCELS L 346 and MCELS L 390, are 
confirmed eclipsing binaries. MCELS L 346 contains a 
nearly equal mass binary, and is modeled with two ioniz- 
ing stars. The rest of our sample show no obvious signs 
of OB companions, such as composite spectral types or 
large discrepancies between Q Q jj q and <2o,Sed, and as 



discussed in S2.2 the binary status will not affect our 
results. 

We evaluate the stellar atmosphere models by match- 
ing the emission-line spectra predicted by CLOUDY 
photoionization simulations with that observed in the 
nebulae. Our nebular simulations show that uniform 
Stromgren spheres cannot reproduce the observations. 
In these cases, [Sn], [On] and [Nil] are under-predicted, 
which suggests that the ionization parameter is too high 
in these simulations. We explore the effects of chang- 
ing i?i nncr and Z, but find that neither parameter can 
resolve this issue. Additionally, we find that [Sn] and 
[Oi] depend strongly on the nebular parameters, but only 
weakly on the particular SED used in the simulation. In 
contrast to the uniform Stromgren spheres, simulations 
that assume a clumpy medium reproduce the observed 
data. The clumpy medium is described by specifying a 
filling factor of dense gas and assuming that the remain- 
ing volume is a vacuum. In general, we find that a filling 
factor of ~ 0.10 produces simulations that agree with 
observations. 

We compare the rate of ionizing photons predicted by 
the best fitting atmosphere, assuming a constant log(g), 
with the rate calculated from L(Ha). Within the errors 
due to uncertainty in log(g), which can be as high as 
a factor of two, we find that the two values are consis- 
tent. However, we do find a systematic offset between 
the rate of ionizing photons predicted by different atmo- 
sphere models. TLUSTY produces the highest Q , while 
WM-basic and SNC02 have slightly lower rates. In con- 
trast, the Qq from CoStar is a factor of two lower than the 
other atmosphere models. The hardness of the CoStar 
SED can explain this offset. We select the best-fitting at- 
mosphere models by matching the ionization balance in 
the nebula. Since the CoStar models have harder SEDs, 
the best-fitting CoStar atmospheres are a few thousand 
degrees cooler and produce fewer ionizing photons as a 
result. 

We evaluate the stellar atmosphere models by compar- 
ing the predicted to observed emission line ratios. For 
emission lines with IP < 36 eV, we find that simula- 
tions ionized by WM-basic and SNC02 atmospheres, re- 
produce the observations. Simulations using TLUSTY 
are also consistent, with the exception of [On] A3727. 
For this line, TLUSTY over-predicts the amount of [On] 
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Figure 12. A comparison of the T e g- from this paper to different SpT-T c ff calibrations. Symbols represent temperatures derived from 
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A3727 in the nebula. Surprisingly, despite having the 
most approximate treatment of line blanketing, simula- 
tions using CoStar reproduce the emission lines up to [O 
III], at 35.1 eV, but they do so at lower T c g than the 
other model atmospheres. 

At energies > 36 eV, the predictions from different 
atmosphere models diverge from each other and the ob- 
servations. We find significant scatter between the pre- 
dictions and observations of [Nem] A3869 from all at- 
mosphere models. CoStar, the hardest atmosphere over- 
predicts the [Nem] line by more than a factor of two in 
most of our objects. In contrast, the other atmosphere 
models range from over- to under-predicting [Nem] by 
70%. Simulations ionized by TLUSTY, with its softer 
SED and plane-parallel geometry, systematically produce 
less [Nem] than simulations using WM-basic and SNC02. 

Finally, we compare the SpT-T e g- calibrations in the 
literature with the best fitting T e g from the simulations 
in this study. Our results for stars with S pT earlier than 
08 V fall along the LMC calibration of |Massey et al. 
(2005), while stars later than 08 V seem to be hotter 
than the calibrations predict. The SpT-T c ff calibrations 
from the literature are based on photospheric lines, while 
the effective temperatures obtained from this work de- 
pend the ionization balance of the nebulae. The com- 
parison of these two methods has the potential to reflect 
differences between the outer atmospheres of the mod- 
els and the actual ionizing stars. However, there are a 
variety of uncertainties and our sample is small, which 
limits its usefulness in this regard. Further study with a 
larger sample and detailed stellar modeling is needed to 
determine if the deviations at late SpT are significant. 

This work uses single-star HII regions to test the pre- 
dictions of stellar atmosphere models. By using single- 
star Stromgren spheres, we dramatically reduce the free 
parameters involved in evaluating atmosphere models 
with observed nebulae. Future directions that could im- 
prove on the groundwork presented here include leverag- 
ing a larger range in T c g by including stars earlier than 
05.5 V, using a detailed fit to high resolution stellar spec- 
tra to generate the input model SEDs, and evaluating 
model atmospheres at different metallicities. 
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